SPECTRAL ZETA FUNCTIONS OF FRACTALS AND 
THE COMPLEX DYNAMICS OF POLYNOMIALS 



ALEXANDER TEPLYAEV 



Abstract. We obtain formulas for the spectral zeta function of 
the Laplacian on symmetric finitely ramified fractals, such as the 
Sierpihski gasket, and a fractal Laplacian on the interval. These 
formulas contain a new type of zeta function associated with a 
polynomial (rational functions also can appear in this context). It 
is proved that this zeta function has a meromorphic continuation 
to a half plain with poles contained in an arithmetic progression. 
It is shown as an example that the Riemann zeta function is the 
zeta functions of a quadratic polynomial, which is associated with 
the Laplacian on an interval. The spectral zeta function of the 
Sierpihski gasket is a product of the zeta function of a polynomial 
and a geometric part; the poles of the former are canceled by the 
zeros of the latter. A similar product structure was discovered by 
M.L. Lapidus for self-similar fractal strings. 
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1. Introduction 

In this paper we obtain formulas for the spectral zeta function of the 
Laplacian on symmetric finitely ramified (or p.c.f., see j2UI22]) fractals, 
such as the Sierpihski gasket and an interval with a fractal measure. 
These formulas involve a new type of zeta function associated with a 
polynomial, and we show that this zeta function has a meromorphic 
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continuation to a strip beyond the half plane where it is analytic. In the 
case of the Sierpihski gasket the spectral zeta function has a product 
structure, with one term which can be described as a "geometric" part, 
and the other term which is the zeta function of a certain quadratic 
polynomial. Both of these features, the existence of a meromorphic 
continuation and a product formula, are reminiscent of the theory of 
one dimensional fractal strings first defined and studied by M. Lapidus 
in [27]. We present basic definitions and results in Section |21 and the 
reader can find more details and references in (2E1 I2H1 EH ED EH EHj- 
In particular, equations (5.2)-(5.3) in [27| (see also section 1.3 of j2H]) 
contain the product formula 

where ((s) is the Riemann zeta function, and 

is the geometric zeta function of a fractal string £, a disjoint collec- 
tion of intervals of lengths The spectral zeta function Cl( s ) °f a 
normalized Dirichlet Laplacian L on L is defined as 

c^) = E A 7 s/2 > 

where Xj are the nonzero eigenvalues of L. 

The first step in our work is to define and study a new type of zeta 
function associated with a polynomial (see Section El)- Suppose that 
R(z) is a polynomial of degree N such that R(0) — 0, c = R'(0) > 1, 
and the Julia set of R{z) lies in IR + . Then we define the zeta function 
of the polynomial R(z) as 



Cr,J s ) = to* E ( C ^)~ S/2 (1-3) 

for 



^e-R- n {z } 



Re(3 ) > dR = (1 .4) 
logc 

where i?~™{2o} means the set of N n preimages of zq under the nth 
composition power of the polynomial R. Moreover, we prove that there 
is an increasing unbounded sequence of positive real numbers Xj such 
that ( R ZQ can be represented by the Dirichlet series 

oo 

W*) = 5>7 S/2 (i-5) 

3=1 
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if fll.4j) holds. We give examples where Xj are actually eigenvalues of 
a Laplacian, but it is possible to define Xj in terms of R and Zq only, 
without reference to any Laplacian. 

Although the definition above make sense under less restrictive con- 
ditions on R(z), they are automatically satisfied if the polynomial is 
associated with a self-similar Laplacian (see (3312111121), as it will be 
the case in the applications we consider. 

Our main result is that ( Rzo {s) has a meromorphic continuation of 
the form 

where £ R (s) is analytic in C and i> Rzo (s) is analytic for Re(s) > 0. If 
3r is totally disconnected, then there exists e > such that i> Rzo (s) 
is analytic for Re(s) > —e. The set of poles of ( Rzo (s) is contained in 
|2]ogA^m7r . n ^1 an j there always is a pole at d R . 

Independently of our work this result was recently (after the resent 
paper was completed) improved in [T3| based on the ideas of an earlier 
paper by Peter Grabner. In particular, it is shown in ^3] that 
Cr 20 actually has a meromorphic continuation to the entire complex 
plain for any polynomial R satisfying the same conditions as in our 
paper. Moreover, [T3| contains a detailed analysis of ( Rzo in the case 
of the Sierpihski gasket. However, an advantage of our method is that 
it is also applicable in the case when R is a rational function satisfying 
appropriate assumptions. Therfore it is likely that our results can be 
extended to the class of fractals considered in 

One of the reasons this new class of zeta functions is interesting is 
that the Riemann zeta function ((s) is the zeta function of a quadratic 
polynomial (Theorem |3J). This result is obtained by comparing two 
representations of the Riemann zeta function: as the Dirichlet series 
and as the spectral zeta function of the standard Laplacian on an inter- 
val. Then the interval can be considered as a self-similar fractal, and 
the spectrum of the Laplacian can be described in terms of the qua- 
dratic polynomial. This construction is not unique in the sense that 
the Riemann zeta function can be represented as the zeta function of a 
polynomial of any degree (see Theorem 01 for an example). In Section |S] 
this construction is generalized for fractal Laplacians on the interval, 
that is, Laplacians associated with a fractal self-similar measure. 

In Section |U] we obtain the following formula for the spectral zeta 
function of the Laplacian A M on the Sierpihski gasket 

Ca m ( s ) = (wiz3 + H^fc) +Ki?„|(' s ) (5573=3 ~ f^=i) ( L7 ) 
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where R(z) = z(5 — Az). 

Combining (J1.6)) and (|1.7|) we obtain that the spectral zeta function 
of the Sierpihski gasket has a form which resembles the spectral zeta 
function of a self-similar fractal string (jl.lj) . For example, for the Can- 
tor self-similar fractal string L, that is the complement of the middle 
third Cantor set in [0, 1], we have 

Cl(*) = C(s)/(3 s -2). 

Thus the spectral zeta function of the Sierpihski gasket has the same 
structure as the spectral zeta function of a disjoint collection of "frac- 
tal intervals". However, this analogy is somewhat superficial because 
such "fractal intervals" so far have not been constructed or defined. 
This means that the spectrum has a product structure which does not 
manifest itself in the structure of the underling space. 

The poles of a spectral zeta function are called the complex spectral 
dimensions (see [271I2E1I2II)- They are 1 and { ^f+f^ : n E Z} for the 
Cantor self-similar fractal string, and {gf , log ^ n7T : n G Z} for the 
Sierpihski gasket, see Figures l2~2l and An interesting feature of the 
product formula (jl.7|) is that, in the case of the Sierpihski gasket, the 
poles of the "polynomial" part of ( A (s) are canceled by the zeros of 
the "geometric" part of Ca m ( s )- We conjecture that such a cancelation 
takes place for other symmetric fractals, and can be explained by the 
interplay between geometric and analytic properties of the fractals. 

The main results of this paper were announced, without proofs, in 
|¥T] with a few typos, which we correct in the present paper. 

Acknowledgment. The author is deeply grateful to M. Lapidus for men- 
toring, encouragement and support, and to R. Strichartz and J. Kigami 
for many helpful discussions. The author very much appreciates inter- 
esting and useful remarks and references by P. Grabner. The author 
thanks the anonymous referee for many helpful suggestions. 

2. The spectral zeta function of a self-similar string 

In this expository section we follow the work of Lapidus et al. (see 
EH |2H1 123 EOl ED E2 and the references therein). 

Definition 2.1. If A is a nonnegative self-adjoint operator with dis- 
crete spectrum and positive eigenvalues {A.,}, counted with multiplici- 
ties, then its spectral zeta function is 



a(*) = £A7* 
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Note that zero is not included in the sum, even if it is an eigenvalue of 
A. 



Example 2.2. If A=— -I3 is a Neumann or Dirichlet Laplacian on [0, 1] 
then 



00 /.2VS/2 



i=i 

where £(s) is the Riemann zeta function. The series converges for 
Re(s) > 1, and has a meromorphic continuation to C with one pole at 
s = 1. 

Definition 2.3. The poles of Ca(s) are called the complex spectral 
dimensions of A. 

Definition 2.4. If £ is a fractal string, that is, a disjoint collection of 
intervals of lengths 4, then its geometric zeta function is 

Cc(«) = £ Z *- 

The next theorem is due to M.L. Lapidus [23 EHl EH1 EDI EH E21 ESJ • 

Theorem. If A = —-^2 is a Neumann or Dirichlet Laplacian on £ 
then 

Us) = ^C(s)Us). 

Example 2.5 (a self-similar fractal string). Let r > 1 and lj = cr~ k 
for N h intervals, where iV is a positive integer, k = 0, 1, 2, ... . Then 

= , c * > Ca(*) = C00: 



If < r then the total length of £ is finite and the spectral dimension 
is 1, otherwise it is 7^7- • The complex spectral dimensions are 1 and 

L log r J 

For example, for the Cantor self-similar fractal string, that is the 
complement of the middle third Cantor set in [0,1], we have c = |, 
r = 3,N = 2. 



I 2 

3 3 



Figure 2.1. Cantor self-similar fractal string. 
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Figure 2.2. Complex spectral dimensions of the Cantor 
self-similar fractal string. 

3. The zeta function of a polynomial 

We suppose that R(z) is a polynomial of degree N such that R(0) = 
and c = R'(0) > 1. We also assume that the Julia set 3r of R(z) 
satisfies Sr C K+. These assumptions imply immediately that the 
coefficients of R are real. Since zero is a repulsive fixed point, its 
preimages are dense in J R by the classical paper [TT]. Thus, Sr C M+ 
if and only if there is a real interval [0, a] C IR+ such that -R _1 [0, a] C 
[0,a]. 

By [I^, the Julia set in our situation is an interval or a Cantor set of 
Lebesgue measure zero. We denote the convex hull of Sr by I, which is 
an interval. By the same paper, if 3r, is an interval, then R~ n (I) = I. 
Otherwise R~ n (I) QI / consists of N n disjoint intervals. In both cases 
R" n -\I) C R-' n {I) and 

00 

J R =f]R- n (I), 

n=l 

where R~ n (-) denotes the preimage under the nth composition (itera- 
tion) power of the polynomial R(z). Also one has a unique decomposi- 
tion / = Uj=i h-> wnere Ij are intervals with disjoint interior, and R(z) 
is monotone on each Ij. 

In what follows we denote I = f]f =1 R{Ij)- Then we have -R -1 (/ ) Q 
I C / . 

Definition 3.1. The zeta function of the polynomial R(z) is 
<R,J S ) = lim E ( C ^)" S/2 

zeR- n {z }\0 
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where z G J , and Re(s) > cZr. Here we define 

21ogiV . , 

d R = — -2—. 3.2 
logc 

Theorem 1. For i?e(s) > i/je limit h3.1\) exists and is an analytic 

function of s. Moreover, there is a nondecreasing unbounded sequence 
of positive real numbers Xj such that 

oo 

w*) = E A 7 s/2 ( 3 - 3 ) 

3=1 

where Xj are defined in terms of R and zq. 

Proof. By the monotonicity, the restriction of the polynomial R(z) to 
each Ij has an inverse, which we denote by Rj . We assume that the 
intervals Ij are numbered consecutively starting from the origin. In 
particular, R^ n (-) is the nth composition power of the branch of the 
inverse function of R(z) that passes through the origin. 

We denote by d+ the part of 3r that is not in R± (I), that is 

TV N 

d + = 3 R 0([jR- k 1 (I)) = [jR k 1 (3 R ). 

k=2 k=2 

We denote by 51 the limit 

%(z) = lim c n R^ n (z). (3.4) 

n^oo 

Since R'(0) = c, this limit exists for any z G Iq and satisfies 5l(z) = 
c'R(R^ 1 {z)) and is an increasing continuous function. The function 31 
was first defined and studied by P. Fatou in [T3] . Its analytic properties 
are well studied in complex dynamics, but we need only the elementary 
facts that the limit (j3.4j) exists and is an increasing continuous function 
for real z G Iq. To show this, notice that R'(z) > and R"(z) < 
in the interior of I\ by 11 . Then it is easy to see that z n = R± n (z) 
converges to zero exponentially fast. Hence 

cz n+1 = cR^(z n ) = z n + O(z 2 n ) n __ too = z n (l + 0(z n ) n _ ¥ 

which proves ([3.4)1 . 

Let W n be the set of sequences (words) w — k± . . . k n of length n, 
where hi G {1, . . . , iV}. For a word w — k\ . . . k n G W n we write _R~ n = 
R k ^ . . . R^ . We use the convention that Wq = {0} and R^{z) = z for 

W = G Wq. 
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Let Wj~ be a subset of W„ of words that do not start with 1. Then 



we have 



n 



(c n z)- s/2 



C R ,z ( S ) = lim E 
zeR- n { 

n 

= lim V c- sm/2 V (c n - m R^ n+m RZ m {z )) 



„-, x ^ 2_. ( c ^ W (3.5) 

m=0 ™ew+ 



m=0 w£W+ 



s/2 



for Re(s) > d/j. 

Here and in what follows, in order to simplify notation, we use the 
convention that z~ s l 2 = for all real s if z = 0. 
The last expression in ()3.5|) can be rewritten as 

oo oo 

c«» = E E (^ m (-o))r /2 = E A 7 s/2 ( 3 - 6 ) 

m =°w€W+ i =1 

where the nondecreasing unbounded sequence of positive real numbers 
{Xj}'jL 1 is obtained by reordering the countable locally finite set of 
positive real numbers Um=o LL e w+ "l "^ G^w^O 2 *)))} \{0}, taking the 
multiplicities into account. □ 

Theorem 2. The zeta function Cr, Zo ( s ) °f ^ e polynomial R(z) has a 
meromorphic continuation of the form 

W (^)= 1 _^ c l /2 +^ ( g ) 

where £ R (s) is analytic in C and i> Rzo (s) is analytic for Re(s) > and 
is defined by (|3.7j) . T/ie se£ of poles of ( Rzq (s) is contained in 



and i/iere always is a pole at d R . 

U3r is totally disconnected, then there exists e > such that ip Rz (s) 
is analytic for Re(s) > —e. If, in addition, r = max < 1 then 

z£3r 

we can take 

2 logr 
log c 

Independently of our work this result was recently improved in [T3j . 
In particular, it is shown in ^3] that ( Rzo has a meromorphic con- 
tinuation to the entire complex plain for any polynomial R satisfying 
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the same conditions as in our paper. However, an advantage of our 
method is that it is also applicable in the case when R is a rational 
function satisfying appropriate assumptions. Therfore it is likely that 
our results can be extended to the class of fractals considered in 



Proof of Theorem® To the last expression in (|3.5|) we add and sub- 
tract 



where 



Us) 



1 - Nc~ s / 2 



{R{z))~ s/2 q!k{z). 



(3.7) 



Here we use the fact (see, for instance |H1 El HI] ) that on 3r there 
is a unique balanced invariant probability measure k which satisfies 



for any sequence j±, j m . 

Note that £ R {s) is an analytic function since 3+ is separated from 
zero, and that l _ N 1 c - s / 2 — Em=o J\f m c~ sm / 2 . Then we only have to show 
that 



-sm/2 



m=0 
oo 



(fwc 1 ^)) 



-s/2 



N m / (a(z))-'/ 2 dK(z) 



m=0 



s/2 



s/2 



dniz) 



Rw m {3+} 



is analytic when Re(s) > 0. 

Thus we have to show that the series converges uniformly when s is 
in a compact subset. To this end we can estimate the absolute value 
of the expression in square brackets by 

M max \z-R- m {z )\^MX{R- m {I Q }) 



where A is the Lebesgue measure on K, 



M = max 

zei+ 



d 

dz 



(*(*)) 



-s/2 



and J + = UweWi+^Uo}- 
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Then, since K(R- m {3}) = N~ 



-in 



j |(^ m (^y 



3,2 - {%{z))- s/2 



cIk(z) ^ 



MN- m X(R- m+1 (I)). 

In this estimate we have to use the exponent — m + 1 instead of m 
because we assumed that z G Iq, not necessarily z G I. 

Thus we have the required convergence for Re(s) > 0. 

If 3r is totally disconnected, then by pages 122-124 in [TT], there is 
r < 1 such that 

X(R- n (I)) < const -r n . 

This estimate gives the required convergence for Re(s) > — e = 2 lo ° g c r . 
Finally, note that r max □ 

z£3r 

In the applications of this theorem considered in this paper, the 
assumptions on the polynomial R(z) given above are always satisfied 
because the Julia set of R(z) is contained in the spectrum of a nonneg- 
ative self-adjoint operator (see [SU E2 ES] ) • Moreover, 

w*)=E A r /2 ( 3 - g ) 

where {Xj} are eigenvalues of a Laplacian, and so the spectral zeta func- 
tion has a representation as a zeta function of a polynomial. In these 
cases the parameter z usually depend on the boundary conditions. 
However, even if the polynomial R(z) is not related to any Laplacian, 
the representation (J3.8|) still holds with {Xj} defined in terms of the 
limit points of [j n>0 c n R~ n (z ) as in Theorem [T] 

Theorem 3. The Riemann zeta function ((s) has a representation 

where C = ^/2tt and ( R0 (s) is the zeta function of the polynomial 
R(z) = 2z{2-z). 

Proof. For the Neumann or Dirichlet Laplacian A on I = [0, 1] we have 



CO 



o=i i=i 

Let 
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be the probabilistic Laplacian on {0, . . . , ^±1, . . . , 1}, which is the 
generator of the simple nearest neighbor random walk. Then 

A/(x) = -f"{x) = 2 lim 4 n A n f{x) 

for any twice continuously differentiable function /. The theorem then 
follows from the definition of Cro( s ) ano - Lemma fa .21 since 

CaOO = l^^UA 8 ) = \i s/2 U^)- 

Note that the preimages of zero are {0,2} = cr(A ), and so the spec- 
trum of A n as a set coincides with i? _n_1 (0), which explains the extra 
factor 4 S//2 above. Also note that 1 is the preimage of 2 of multiplicity 
two, but 1 as a simple eigenvalue of A n , which explains the extra factor 
| above. □ 

Lemma 3.2. If z ^ 1 then R(z) = 2z(2 — z) G a(A n ) if and only 
if z G cr(A„ + i), with the same multiplicities. Moreover, the Neumann 
discrete Laplacians have simple spectrum with cr(A ) = {0, 2} and 1 G 
a(A n ) for all n > 0. 

Proof. This lemma can be proved by an elementary computation, see 
[311 Example 2.4 and Theorem 5.8]. □ 

4. Mellin transform 

In this expository section we contrast the zeta function of a polyno- 
mial with the Mellin transform of an invariant measure related to the 
same polynomial. We do not use the correspondence between these 
object in our work, but the reader can find some related formulas in 

Definition 4.1. If k is the unique balanced invariant measure on the 
Julia set 3r of R(z) then its Mellin transform is 

MJs) = lim N~ n V z- s/2 = [ z- s/2 dK{z). 



z&R-"{z } 



3r Ai A 2 A 3 

H - - - - e e s 

Figure 4.1. A pictorial illustration of the Julia set 3r 
of R(z) and the sequence {Xj} in Theorem [T] 



12 



ALEXANDER TEPLYAEV 



JLr 



dp. 



o o 

o o 

o o 

o o 

o o 

d R + 1 d R + 2 



o o 

o o 

o o 

o o 

o o 



Figure 4.2. A pictorial illustration of the poles of the 
zeta function of a polynomial in Theorem |2] (left), and of 
the poles of its Mellin transform (right). 



Theorem (7\). The Mellin transform M R converges for Re(s) < oIr 
and has a meromorphic continuation to C with the set of poles contained 
in 

' 2 log N + Ain-n 

\- m : n E 1*, m G N 

logc 

Every real pole has a positive residue and thus is not degenerate. 

Note that by the Theorem of Bessis, Geronimo and Moussa the 
Mellin transform has infinitely many nondegenerate positive real poles, 
while by Theorem|2]the zeta function of a polynomial has only one pos- 
itive real pole (see Figure | 



5. Fractal Laplacian on an interval and its zeta function 

We begin with defining the usual energy form of the unit interval 
by introducing a self-similar structure. In this way the interval can be 
considered as a p.c.f. fractal (see [21 EH [2H 1221)- In particular, in these 
papers the reader can find proofs of the propositions given below, as 
well as the general theory of Dirichlet forms on fractals and further 
references. 

We use three contractions because it will allow us to easily modify 
the construction in order to obtain a fractal Laplacian on the unit 
interval. One can make the same construction using any number of 
contractions. Note that in Theorem El we essentially used a self-similar 
structure on the unit interval with two contractions. 
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Remark 5.1. We emphasize that in this and the next section we con- 
sider two different iteration function systems which are not to be con- 
fused. We start with an iteration function system of linear contractions 
to define a self-similar set, which is the set where the Laplacian is de- 
fined (and the state space of the corresponding diffusion process). This 
self-similar set is often a fractal, but in this section it is an interval 
with a self-similar structure. Then, based on the "spectral decima- 
tion" or, more appropriately called, "spectral self-similarity" (see 
we describe the spectrum of the Laplacian using inverse iterations of a 
polynomial (a rational function in more general situation of |34j). To 
this latter iteration function system we apply the results of Section El 

We define three contractions Fi, F2, F3 : M 1 — > M. 1 with fixed points 
xj = 0, |, 1 by 

Then the interval I=[0, 1] is a unique compact set such that 

3=1,2,3 

The discrete approximations to I are defined inductively by 

v n = U*i(v^) = {£}L 

j=l,2,3 

where Vq = dl = {0, 1} is the boundary of /. For x, y £ V n we write 
y ~ x if \x — y\ = 3~ n . 

Definition 5.2. The discrete Dirichlet (energy) form on V n is 
UfJ) =3"£(/(y)-/(x)) 2 

where the sum is taken over all unordered pairs x, y £ V n such that 
y ~ x. The Dirichlet (energy) form on / is 

£(/,/)= lim £„(/,/) 

if this limit exist and is finite. By the next proposition this limit is 
increasing and so it always exists as a nonnegative real number or +00. 

A function h is harmonic if it minimizes the energy given the bound- 
ary values. 

Proposition 5.3. We have that £ n+ i(/, /) ^ £ n (/, /) for any function 
f , and 

£ n +i(h, h) = E n (h, h) = E(h, h) 
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for a harmonic h. A function h is harmonic if and only if it is linear. 
If f is continuously differentiable then 

£(/,/)= f\f{x)\ 2 dx. 
Jo 

The Dirichlet (energy) form £ on I is self-similar in the sense that 
£(/,/) = 3 £ £(/oF,,/oF,). 

j=l,2,3 

Definition 5.4. The discrete Laplacians on V n are 



yev n 



and the (renormalized) Laplacian on / is 



A/(x) = lim TA n f(x) = \f"{x) 



for any twice differentiable function. 

Obviously, for any twice differentiable function this Laplacian satis- 
fies the Gauss-Green (integration by parts) formula 



£(/,/) = -2 f fAfdx + ff 
Jo 



and has the spectral asymptotics 

hm — —- = - 

A->oo A ds / 2 7T 

where p(A) is the eigenvalue counting function p(A) = #{j : Aj < A} 
of the Dirichlet or Neumann Laplacian A and d s = 1 is its spectral 
dimension. 

In what follows we modify this construction to obtain a one param- 
eter family of fractal Laplacians on the unit interval. We introduce a 
parameter < p < 1 and define q — 1—p, which later will be shown to 
have a meaning of transitional probabilities of a random walk. Then 
we define contraction factors (or "resistance weights" ) 

p + pq 2q — pq 

ri = r 3 = — and r 2 = — — , 

2 + pq 2 + pq 

and "measure weights" 

q p 

mi = m 3 = and m 2 = . 

1 + q 1 + q 



SPECTRAL ZETA FUNCTIONS OF FRACTALS AND POLYNOMIALS 15 

Note that, up to a constant factor, the resistance weights are reciprocals 
of the measure weights, and 

mi + m 2 + m 3 = r 1 + r 2 + r 3 = 1. 

By their definition, constants resistance weights and measure weights 
rrii depend on the choice of p. 

We define three contractions F Pt i, F Pj2 , F Pi3 : IR 1 — > M 1 with fixed 
points Xj — 0, |, 1 by 

F pA x ) = r j x + ( 1 ~ r j) x j- 
Then the interval I—[0, 1] is a unique compact set such that 

3=1,2,3 

The discrete approximations to / are defined inductively by 

3=1,2,3 

where, as before, V = dl = {0, 1} is the boundary of /. 

Definition 5.5. The discrete Dirichlet (energy) form on V n is defined 
inductively 

£»»(/,/)= E :(. /•-/•:,)• 

3=1,2,3 

with £ (/,/) = (/(l) - /(0)) as before. The Dirichlet ( energy) form 
on / is 

E p (f, f) = lim E Ptn (f, f) 

if this limit exist and is finite. Note that this is the same definition 
as the Dirichlet form defined before if and only if p = ~. By the 
next proposition this limit is increasing and so it always exists as a 
nonnegative real number or +oo. 

The next two propositions follow, for example, from the general re- 
sults in |21 EH El 122] • Again, a function h is harmonic if it minimizes 
the energy given the three boundary values. 

Proposition 5.6. We have that E Pin+ i(f, f) ^ £ p , n (f, f) for any func- 
tion f , and 

£ p ,n+i{h, h) = E Pt n(h, h) = E p (h, h) 

for a harmonic h. 
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Proposition 5.7. The Dirichlet (energy) form £ p on I is local and 
regular, and is self-similar in the sense that 

£*(/./)= E -^,i/o/-;, r /o/-;,. ; ). 

J'=l,2,3 

T7ie domains of£ p and of the corresponding Laplacian A M are contained 
and dense in the space of continuous functions. 

There exists a fi-Laplacian A M , that is a densely defined self-adjoint 
operator that satisfies the following Gauss-Green ( integration by parts ) 
formula 

e P (/,/) = -c f fAjdn + ff'H. 

Jo 

where /i is a unique probability self-similar measure with weights mi, 
7712, ^3, that is 

/i = rrij fioFpj. 

5=1,2,3 

A,J(a;) = lim(l+A)"A Pin /(x) 
where the discrete Laplacians 

{pf{yh-x) + qf(yk+i) - f(yk) 
or 
qf(vk-i) +pf(yk+i) - fivk) 

are defined as the generators of the nearest neighbor random walks on 
V Pt n with transitional probabilities p and q assigned according to the 
weights of the corresponding intervals. 

Note that by definition p = — ^ — , q = — ^ — . The transitional 

J 1 mi+m2 ' ^ mi+m2 

probabilities p and q can be assigned inductively as shown in Figure I5~T1 
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FIGURE 5.1. Random walks corresponding to the dis- 
crete Laplacians A p n . 
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Note that the previous construction of the standard Laplacian and 
the Dirichlet form on / = [0, 1] corresponds to p = |. If p ^ | then 
we can make a change of variable on the unit interval / so that either 
the Dirichlet form becomes the standard one, or the measure becomes 
the Lebesgue measure, but not both. Thus for different values of p the 
/i-Laplacians are different even up to a change of variable. 

By the result of Kigami and Lapidus [23], the Dirichlet or Neumann 
Laplacian A M has the spectral asymptotics 

Pp(x) PvW 

< hm inf , , j- < hm sup - , < oo 

where p p (X) is the eigenvalue counting function of A M , and its spectral 
dimension is 

log9 , 

" log " ' 

where the inequality is strict if and only if p ^ q. 
In what follows we denote 

R p (z) = z{z 2 + 3z + 2 + pq)/pq. 

We use notation cr(A) for the spectrum of a linear operator A. 

Lemma 5.8. If z ^ —l±p then R(z) G a(A pn ) if and only if 
z G cr(A pra+1 ) ; with the same multiplicities. Moreover, the Neumann 
discrete Laplacians have simple spectrum with cr(A ) = {0, —2} and 

n-l 

a(A p , n ) = {0,-2}|J VH" 1 ^} 

m=0 

for all n > 0. In particular, for all n > we have 0, — 1 ± q, —2 G 
<j(Ap jn ). v4/so, for all n > we /iai>e — 1 ± p G cr(A Pin ) and on/y if 
p = q. 



Proof. In this case, according to jJUJ Lemma 3.4], [3U (3.2)] or 
we have that R„(z) = ^ X . . where o9 and o?i solve the matrix equation 

S-zIo- X(Q - zhy'X = vo(z)H - ifxWo- 
with S = X = X = Iq = ii, the 2x2 identity matrices, 

— 1 p 



Q 

and 



P 
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Solving this we obtain 

pq 



z 2 + 2z + 1 - p 2 
and 

( \ _ z(z 2 + 3z + 2 + pq) 
MZ) ~ z 2 + 2z + l-p 2 ■ 
Then we use the abstract spectral self-similarity result (see |4U[ |33J 021 
HH]) that the spectrum of A p n+1 can be obtain from the spectrum of 
by the inverse of R(z). Note that and 2 are fixed points of R(z). 
The preimages of are 0, —1 — p and — 1 — q. The preimages of 2 are 
2, — 1 +p and — 1 + q. If p ^ q then — 1 ±p are not eigenvalues because 
they are poles of y>„(z) (see (331110111211331) □ 

Theorem 4. If p ^ \ then d s = oIr < 1 and t/ie spectral zeta function 
of is 

ca» = iz |^ (w«)+Ch, w oo) 

where c = 1+-^ owe? 21,22 = 1 ± g. Moreover, Ca m ( s ) ^ as a mesomor- 
phic continuation to Re(s) > —e for some e p > 0, and liminf e„ ^ 2, 

p— >o 

lim inf £„ > 2 . 

Proof. This follows from Lemma 15. 8^ Theorem 121 and the definition of 
the zeta function of a polynomial. The last statement follows from the 
fact that as p —>■ or p — > 1 we have limmin \R'(z)\ = oo. □ 

Theorem 5. The Riemann zeta function ((s) has a representation 

COO = \cx R)0 (s) 

where C = % and Cro( s ) ^ s ^ e ze ^ a function of the polynomial 

R(z) = z(Az 2 - 12,2 + 9). 

Proof. The proof follows from Lemma 15 .tt\ Theorem 121 and the defini- 
tion of the zeta function of a polynomial, analogously to the proof of 
Theorem El and of Theorem 0] with p = q. □ 

Remark 5.9. In Figure I5~2*l we give a sketch that describes the complex 
dynamics of the family of cubic polynomials associated with the fractal 
Laplacians on the interval. Note that the Julia set is real, and so what 
we sketch is actually the real dynamics of these polynomials. The 
curved dotted line corresponds to the case when p = \ and the Julia 
set of R(z) = Ri(z) is the interval [—2,0]. For any other value of p 



2 
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-2 



Figure 5.2. Sketch of the cubic polynomial R p (z) as- 
sociated with the fractal Laplacians on the interval. 

(0 < p < 1), the graph of the polynomial R p (z) behaves like the shown 
solid curved line. It is easy to see that then the Julia set of R p (z) is a 
Cantor set of Lebesgue measure zero. Note that the change p h- > 1 — p 
does not change the polynomial R p (z), although the Laplacians A M are 
different. 

6. Spectral zeta function of the Sierpinski gasket 

We begin this section with describing the analysis on the Sierpinski 
gasket in the way first developed by J. Kigami (201 • We recall only 
basic facts here, and in the papers 011211011121111111011131211122112211221 

EUl EZl ED] the reader can find proofs of the propositions given below, as 
well as the general theory of Dirichlet forms and diffusions on fractals, 
many physical applications and further references. 

If we fix three contractions Fj(x) = \{x + Vj), where Vj are vertices 
of a triangle in M 2 , then the Sierpinski gasket is a unique compact set 
S such that 



S 



F 1 (S)UF 2 (S)UF 3 (S). 




Figure 6.1. Sierpinski gasket. 
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One can see that the analysis on the Sierpihski gasket is analogous in 
many respects to the analysis of the fractal Laplacian on an interval we 
described in Sectional This analysis does not depend on the particular 
embedding of the Sierpihski gasket in M 2 , and so a particular choice of 
the points Vj does not matter as long as they are not collinear. 

Definition 6.1. For each n ^ we define V n by 

V n+1 = F 1 (V n )\jF 2 (V n )\jF 3 (V n ) 

where V = {t>i, ^2,^3} = OS is the set of vertices of S. 




Figure 6.2. Discrete approximations to the Sierpihski gasket. 

On each V n there is a notion of neighboring points, giving V n a nat- 
ural graph structure, which can be defined inductively. In Figure 16.21 
the neighboring points are connected by line segments, and in what 
follows are denoted by the symbol u x ~ y" . 

Definition 6.2. The discrete Dirichlet (energy) form on V n is 

£»(/./) = (!)" E (f(y)-f( x » 2 

y~x 

and the Dirichlet (energy) form on I is 

£(/,/)= lim £„(/,/) 

n— >oo 

if this limit exist and is finite. By the next proposition this limit is 
increasing and so it always exists as a nonnegative real number or +00. 

A function h is harmonic if it minimizes the energy given the three 
boundary values. 

Proposition 6.3. We have E n+ i(f,f) ^ £„(/,/) for any function f, 
and 

E n+ i(h, h) = E n (h, h) = E(h, h) 

for a harmonic h. 

Definition 6.4. The discrete Laplacians on V n are 

An/(x) = \ £ f(y)-f(x), xeV n \V 

y&V n 



£(/,/) 
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and the Laplacian on S is 

A,f(x) = lim 5 n A n f(x) 

n— >oo 

if this limit exists and A^f is continuous. 

The next theorem is essentially due to J. Kigami [201 EU 122], al- 
though many relevant results are contained in HT) | IT7 | ITT ^ H3 \ I25j . 

Theorem. & is a local regular Dirichlet form on S which is self-similar 
in the sense that 

£(/./) = ! E ^FjJoFj). 

J=l,2,3 

In particular, the domains of £ and A M are contained and dense in the 
space of continuous functions. The Laplacian satisfies the Gauss-Green 
( integration by parts ) formula: 

-6 j fA fl fd^+J2fiP)df(p) 

where \x is the normalized Hausdorff measure, which is self-similar with 
weights ~, |, | 

J=l,2,3 

anddf(p) is a certain derivative at the boundary. The (Dirichlet, Neu- 
mann ) Laplacian has discrete spectrum and compact resolvent. Green's 
function g(x,y) is continuous. If p(A) is the eigenvalue counting func- 
tion of the Laplacian A^, then 

(J < lim mi - , , n < lim sup - , , n < oo 

where the spectral dimension is 

l<d s = l m<2. 

" log 5 

Next, we present two more theorems illustrating the "fractal" prop- 
erties of the Laplacian A M . 

Theorem (jU The fi-heat kernel of the Laplacian A M (that is, 

the transition density pf (x , y) of the ^-diffusion) has the non-Gaussian 
'fractal" asymptotics 

rf(x, y)~C t~P exp { -C r a d s (x, yf] 
where a= , ffi 9 > /3=j^| and 7 = , 9 . 

log 5— log 2' " log 5 ' log 5— log 2 

Theorem ( jH])- € DomA^ t/ien / 2 ^ DomA M unless f is constant. 



22 



ALEXANDER TEPLYAEV 



Our main result in this section is the following theorem. 

Theorem 6. The zeta function of the Laplacian on the Sierpinski gas- 
ket is given by 

^a m ( s ) = Co |0)-^— • ( i_3.5-s/ 2 + 1-5 3 -/ 2 ) + 



where R{z) = z(5 — 4z). In particular, ( A (s) has a meromorphic 
continuation for Re(s) > —e where e^l. The poles are contained in 

[ 2inix log 9 + 2m7r 

< , : n G Z 

L log 5 log 5 

We have < ^ = d R < I < d s = ^ < 2. 

logo 11 " logo 

Proof. By [T7J HU] , every eigenvalue of has the form 

A = -5 m lim 5 n R~ n (z ) 

where R~ n (zo) is a preimage of zq = 3/4,5/4 under the n-th iteration 
power of the polynomial R(z). The multiplicity of such an eigenvalue 
is of the form C\3 m + C2 where C\ and C2 are computed explicitly 
in [TTl !4Uj . More precisely, by QT). Proposition 3.12], for the discrete 
Laplacian we have that the multiplicity of z E R~ m (|) is 3 " m 2 +3 
for n ^ 1 and ^ m ^ n — 1, and the multiplicity of z G R~ m (|) 

is 3 " ™ 1-1 for n ^ 2 and ^ m ^ n — 2. Then the proof follows 
directly from Theorem El and the definition of the zeta function of a 
polynomial. □ 



Remark 6.5. We have that the zeros of the "geometric part" cancel all 
the poles of the zeta function of the polynomial in the product formula 
(J6.1)) . We conjecture this also for other symmetric p.c.f. fractals [3*4*j . 
except one dimensional. In Figure lfi~3~l the series of the canceled poles 
is denoted by black dots, and the other two series of poles are denoted 
by small circles. 

Remark 6.6. Note that the zeta function of the polynomial in 1)6.1)) 
also describes the distribution of the eigenvalues corresponding to the 
non localized eigenf unctions (see 0110]). 

Remark 6.7. The method of Theorem El applies to a large class of 
fractals (see, for instance, El HOI iH 112 i3 ) , and possibly to |3E|. 
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